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Abstract 

We study the one-pion events produced via neutral current (NC) and charged 
current (CC) interactions by the atmospheric neutrinos. We analyze the ratios 
of these events in the framework of oscillations between three neutrino flavors. 
The ratios of the CC events induced by to that of the NC events and a 
similar ratio defined with help us in distinguishing the different regions of 
the neutrino parameter space. 
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It has been more than a decade since the atmospheric neutrino anomaly has been ob- 
served by the water Cerenkov detectors 1MB and Kamiokande [0,^ . Eventhough conven- 
tional detectors Frejus |^ and NUSEX |^ have not observed this effect, a recent experiment 
SOUDAN-II did indeed see a deficit Recent results from Superkamiokande confirmed 
the earlier results including the zenith angle dependence of the deficit in the multi-GeV data 
1^. The results of all these experiments are presented in the form of the double ratio 

The measured value of R for Kamiokande is O.GOtaoe i 0-05 for sub-GeV data (E < 
1.33 GeV) and 0.57^^? ± 0.07 for the multi-GeV data (E > 1.33 GeV) whereas that of 
Superkamiokande is 0.61 ± 0.06 ± 0.05 for sub-GeV and 0.67 ± 0.06 ± 0.08 for the 
multi-GeV data. The reason for presenting the result in the form of the double ratio is that 
the theoretical calculations of and z/g fluxes are subject to large uncertainties of the order 
of 30% [HH. Hence a comparison of measured neutrino flux with the calculated one (either 
for or for z/g) is not particularly useful. However, the uncertainties in the ratios of the 
calculated fluxes are much smaller (less than 10%). Hence a comparison of the measured 
ratio of flux to Ve flux to the calculated ratio will yield meaningful information on neutrino 
properties. 

Neutrino oscillations provide a natural explanation for both the overall deficit and the 
zenith angle dependence |1T0|, 1 1 1 . Kamiokande analyzed their data in terms of two flavor 



oscillations between z/g and *->■ v^-. In both cases they found allowed regions of 

parameter space, with the mass-squared difference around Am^ ^ 0.02 eV"^ and the mixing 
angle near its maximal value of 7r/4 0. The analysis of Superkamiokande yields roughly 
the same mixing angle but a lower value of Am? ~ 0.001 eV"^ 0)0- Since both ^ 
and ^ Vr oscillations give equally good allowed regions, one must look for ways of 
distinguishing which type of oscillations are the cause of the atmospheric neutrino anomaly. 



Recently Vissani and Smirnov (VS) proposed that the single pion events in the energy 
range 0.5 GeV - 1.5 GeV can be used to distinguish between different types of oscillations. 
Such events induced by neutral current (NC) contain a 7r° whereas the events induced by 
charged current (CC) contain a 7r+ or . The vr^ are detected as sharp rings, whereas 
the 7r° decays into two photons which are detected as two diffuse rings. Kamiokande has 
already observed such vr^'s in their detector by considering all events with two diffuse rings 
and by selecting those events with invariant mass in the range 90 MeV - 180 MeV [jl4 



The CC events due to contain an in addition to a vr^ and those due to z/^ contain 
a fi^. Thus one can distinguish an NC event (two difuse rings with the invariant mass of 
the rings around m^) from a z/g CC event (one diffuse ring due to and one sharp ring 
due to TT^) and from a CC event (two sharp rings from /i^ and vr''^). The two ratios, 
NC events to CC events and NC events to CC events, are free from the uncertainties 
of the theoretical flux calculations. The ratio of the cross sections of the NC to CC events 
in single pion production were measured previously and the uncertainties in them are 
about 15%. Hence a measurement of the above ratios can help in distinguishing the different 
types of oscillations which are currently relevant for atmospheric neutrino anomaly. 

In this paper we analyze the atmospheric neutrino problem in the framework of oscilla- 
tions between three active flavors. We will assume that one of the masses is much greater 
than the other two. We fix the smaller of the mass-squared differences using the solar 



neutrino data [16,17|. Then the larger mass-squared difference and two of the mixing an- 



gles are relevant for the atmospheric neutrino problem |T8[. The double ratio R defined in 
equation(l), along with the zenith angle dependent multi-GeV data, was analyzed in this 
framework previously [p!9| , p0| . Here we analyze the ratio of CC to NC events for one-pion 



production in the framework of three active neutrino oscillations and see how they can dis- 
tinguish between different regions of parameter space. The NC event rate is unaffected by 
the oscillations while the CC event rate is affected by the oscillations and hence depend on 
the neutrino parameters. In section II we give the basic theory for the calculation of the 
one-pion event rates in the three-flavor oscillation scenario. This is followed by the discus- 



sion of the experimentally measurable quantities in section III and a summary is given in 
section IV. 



II. THEORY 

In a three flavor scheme, the weak eigenstates z/q, are related to the mass eigenstates Ui 
through a 3x3 unitary matrix U as 

^a = H^m^i- (2) 
i 

U can be written as 

U = U^\tlj) X UP"""'' X U^\(j)) X U'\uj), (3) 

where, If^^lOij) is the two flavor mixing matrix between the ith and jth mass eigenstates 
with mixing angles 6ij. We assume CP invariance and set t/P^'^*<= = /. For the neutrinos 
that propagate through matter, the CC interaction between z/g and e induces an effective 
mass term for the Ue which is of the form A = 2\/2GFNeE, where Ne is the electron number 
density and E is neutrino energy. Single-pion production maximally occurs for neutrinos 
within the energy range of 0.5 GeV - 1.5 GeV. For this energy region, matter effects in the 
earth can be neglected since A < 10~^ eV"^ and ~ 10^'^ eV"^, as given by Superkamiokande 
analysis [[7|,|r2|. This choice for ^31 is consistent with the mass heirarchy assumption under 
which both the atmospheric and solar neutrino anamalies can be explained. Hence, here we 
have the simple case of vacuum oscillations. The vacuum oscillation probability Pafs for a 
neutrino with a flavor a to oscillate in to a flavor (3 can be given by 



'2.53x^21 



2U^^U^sUp^U,scos (^^:^|^^ + 2U^,U^2Up,U,2Cos (^^:^|^^ . (4) 

Here, E is the energy of the neutrino in GeV, x is the distance travelled by the neutrino 
in kilometers and 6ij is — m| in eV"^. We assume the mass heirarchy (^21 ^ ^32 — 



(531 ) as suggested by the variuos analyses of the solar and atmospheric neutrino problems 



T9|J20|.16-18,2T| and take ^21 — 0, which allows us to put the oscillating term involving 



^21 as unity. So we get 

Pep = {UaiUpi + U^2Up2f + {UM^f + 2{U^iU^i + Ua2Up2)U^3U(33Cos 1^^:!^^ . (5) 

Since we are considering vacuum oscillations with CP invariance, the oscillation probability 
for antineutrinos is the same as that for neutrinos. 

The various one-pion events arising from NC and CC interaction are listed below: 

uiN uiTi'^N, l = e,fi, T, (6) 

viN^l-K^N', / = e,/i, (7) 

VeN e-vr+AT, (8) 

v^N ^ /i-TT+iV, (9) 

where, N and A^' refer to the nucleons. As we are dealing with sub-GeV neutrinos, Ut 
do not excite CC processes. For every neutrino reaction, we consider the corresponding 
anti-neutrino reaction also. 

It is difficult to calculate the zenith-angle dependent fluxes for sub-GeV neutrinos [^,|^ 
and further the scattering angle of the final state charged lepton for neutrinos in the sub- 
GeV range can be as large as 60° and thus directionality is lost in the detection process. 
Moreover, the event rates for the one-pion events that we consider are not very large. Hence 
we will not consider zenith angle dependence and will average over the zenith angle, or 
equivalently, over x, the distance travelled by the neutrinos in its allowed range (20 Km 
-13000 Km). 

The NC rate A^^*-^ (with vr'' as the detectable final state particle) is given by 



rl.bGeV 

AT^^ = / dE 

Jo.5GeV 



cr'^^iE) e^o, (10) 



dE dE 

where, a^^ is the scattering cross section for the reaction of Eq. (^ and e^o is the 7r° 
detection efficiency which is 0.77 for the energy range of interest. A^^*^ being fiavor blind 



remains unaffected by oscillations. The CC rates arising from processes (7-9) can be written 

as 
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a;!(i?)e^oe,-, (13) 

af-{E)e^oe,-. (14) 



In each of the above equations the events due to antineutrinos are added to those due to 
neutrinos. Here, ee+, and e-,,- are the electron, muon and the charged pion detection 
efficiencies respectively, which in principle are energy dependent quantities. However, for 
the limited energy range of sub-GeV neutrinos, they can be taken to be energy independent 
and to be close to unity PJT^. The differential fluxes for electron and muon neutrinos, 
d^e/dE and d(^^/dE, are taken from [p, z/^-nucleon CC cross section is taken from |jl5[ and 
i^e-nucleon CC cross section is obtained by the approximation a^- (E) ~ aj^- {E + 100 MeV) 
0. The NC cross sections for neutrino energy range 0.5 GeV - 1.5 GeV were taken from 
2^ , |25|jr5| . The corresponding rates in the absence of oscillations can be obtained by setting 



the survival probability Paa to unity and transition probabilities Pap {a ^ (3) to zero. 



III. EXPERIMENTALLY MEASURABLE QUANTITIES 

From the processes listed in Eqs. (^ - (|^) we see that a vr*^ is produced both in the 
neutral current process defined in Eq. (^ and in the charged current process defined in 
Eq. (^. The latter process, in general produces three rings (two from the 7r° decay and one 
from the charged lepton) whereas the former produces only two rings. As discussed in the 
introduction, the charged current processes in Eq. (|]) and (^ also produce two-ring events. 
Hence, by choosing only two-ring events we can discriminate against the charged current 
events of the type given in Eq. (|^. However, there is a siginificant probability of about 
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0.17 when the two photons from the 7r° decay cannot be resolved. If the tt" produced in 
Vf.N —>■ eiT^N' is not resolved into two photons, then we see two diffuse rings in this process 
also. Such events can be usually rejected because in general their invariant mass will not 
be near 171.,^. However, if the 7r° produced in z/^A^ fiTc^N' is unresolved, we have a signal 
of one diffuse and one sharp ring and this process forms a background to the reaction in 
Eq. (H). Its effect in all the experimental observables must be properly taken into account. 
An additional process we take into account is the production of two charged pions via the 
neutral current process 

UiN ui n+n-N. (15) 

This process produces two sharp rings and mimics the reaction given in equation (^. We 
take its rate to be about 10% of the rate of reaction in equation (|) |jl3[. 
Let us define the following experimentally measurable quantities: 

1. Ndi:): The number of events with two diffuse rings. 

2. A^^o: The number of events with two diffuse rings whose invariant mass is in the range 
90 MeV - 180 MeV. 

3. Nos'- The number of events with one diffuse and one sharp ring. 

4. Nss'- The number of events with two sharp rings. 

If the detector were ideal and all the particles can be exactly identified, the relations between 
experimentally measured quantities and theoretically calculable quantities would be N^o = 
]\fNC ^ jY^^ _ N^Sl^^ and Nss = ^^^+- The charged particles can be identified with good 
efficiency but the efficiency of identifying vr" events from Nod is only 0.77. There is a 0.17 
probability that the two photons from 7r° decay cannot be resolved and the 7r° appears as 
a single diffuse ring. Thus the CC events with a in the final state lead to contributions 
to Nds and Nss- Hence, the modified relations between the experimental and theoretical 
quantities are 
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N^o = 0.77 AT^^ + 0.17 X OAN^-Zo 



Nss = N^-""^, + N^^^Z- (16) 

In writing the above equations, we added the contribution of u^N n'n^N' to Nds due 
to the 17% tt" misidentification probabihty and the small contribution to A^^^ due to two 
charged pion production. We also took into account the contribution of VeN — > e'Tr^N' to 
A^^o due to the 17% vr^ misidentification probability. Assuming that the directions of the 
electron and the misidentified 7r° are randomly distributed, we estimate that only about 
40% of these events survive the cut on the invariant mass of the two diffuse rings. Hence 
the factor 0.4 in the correction to the expression for A^^o. 

From the experimentally measurable quantities, we define the following three ratios 

= f^. (17) 
R, - ^, (18) 
R, = (19) 

Presently the data on one pion events in the atmospheric neutrino experiments is very 
limited and the errors are quite large. However, we will assume that the large statistics 
of Superkamiokande will enable it to measure all experimental quantities to an accuracy of 
better than 10%. Turning now to theoretical uncertainties, the single largest source of error 



occurs in the ratio 22 



This, when combined with other errors, leads to a minimum theoretical uncertainty of about 
15% in the predictions of -Ri and i?2- Unfortunately, this uncertainty cannot be reduced and 
will be part of all future predictions. Still, even with this handicap, we find that the ratios 
Ri and R2 can be used to distinguish between different types of neutrino oscillations. 



8 



For the atmospheric neutrino oscillations, based on the mass hierarchy assumption and 
using CP invariance, the three relevant parameters are 0, ip and We choose the following 
ranges for the mixing angles: 

• < < 50°. This comes from the analysis of the solar neutrino anomaly. For this 
region of 0, there exist allowed values of uj and ^21 which can account for the various 



solar neutrino data 16,17 . 



30 < ^ < 90°. This comes from our analysis of the zenith angle dependent single ring 
events of Kamiokande data on atmospheric neutrinos. 



For definiteness, we fix ^31 = 0.001 eV^^, which is the central value favored by Su- 
perkamiokande analysis of the single ring atmospheric data []7|,p^. However, the results 



do not depend on ^si significantly. The reason for this is that we are averaging over the 
distance in calculating the number of events and the dependence on ^31 is lost during this 
averaging. We will try to explore how useful the ratios Ri and R2 are in distinguishing 
different regions of the mixing angles. We note that without oscillations (i.e., = and 
■0 = 0), these ratios have the values i?"° = 2.6 and i?2° = 3.4. We plot the values of Ri and 
i?2 as functions of the mixing angle ip for different values of in figures 1 and 2 resectively. 
From these we note that for large values of mixing angles, Ri and R2 differ significantly 
from their no oscillation values. 

Let us consider three particular cases of maximal oscillations: 

1. z/^ ^ z/g oscillations: In this case, if) = 90° and = 45°. Since flux of z/^ is larger than 
that of z/g, the total flux of u^, increases and that of z/^ decreases. Hence Ri increases 
and R2 decreases. For this case, -Ri = 3.8 and R2 = 2.5. Clearly the change in Ri 
and R2 is more than the 15% uncertainty in the theoretical prediction and one can 
distinguish this case from that of no oscillations. 

2. i^fi z/t- oscillations: In this case, ip = 45° and = 0°. Here, the flux of goes down 
and that of u^. is unaffected. Hence Ri does not change and R2 decreases. Numerically 



we have, -Ri = 2.6 and R2 = 1.9. These numbers are significantly different from 
those of the previous case as well as that of no oscillations. Hence the 15% theoretical 
uncertainty is no bar in distinguishing this case. 

3. Maximal three flavor oscillations: Here all the three flavours mix maximally with one 
another. This occurs for ip = 45° and ~ 35°. For this case, Ri = 2.6 and R2 = 2.2. 
Thus the values of -Ri and R2 together allow us distinguish it from the no oscillation 
case as well as the two previous cases. 

Finally, let us consider the ratio R3 which is the sum of Ri and R2 because already one 
can extract a value for it from Kamiokande data. The experimental uncertainty is quite 
large but we will derive the constraints on the neutrino parameter space from the existing 



data. From [|T4[, we estimate R^ = 4.2 ± 0.7. The experimental error is about 16%. To 
this, we add the 15% theoretical error in quadrature. Thus we have R3 = 4.2 ± 1.0. Fig. 3 
shows the plots of R^ as functions of the angle ip for different values of (p. We note that the 
upper limit -R3 < 5.2 imposes the constraint ip < 70°. Hence pure ^ v^, oscillations are 
disfavoured at la level. Figure 4 shows the allowed region in cp — ip plane by the Kamiokande 
constraint on R^. Already, with the very limited data available, the ratios of 1-pion events 
yield useful constraints on the allowed neutrino parameter space. 



IV. SUMMARY AND DISCUSSION 

The atmospheric neutrino experiments have already given us valuable information on 
the neutrino parameters. So far the bulk of the analysis of the experimental data has been 
restricted to the single ring events. With the accumulation of more abundant data in the 
existing and the upcoming detectors, analysis of two-ring events is also likely to be taken up 
in the future. We have used the framework of the three-flavor oscillations, to suggest such 
an analysis. We have shown that independent information on the parameters and ip can 
be obtained. 
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Recently the results of CHOOZ experiment led to strong constraint on the mixing angle 
(f) < 12.5" [p^-|26[]. This in turn implies that the atmospheric neutrino anomaly is almost due 
to z/^ ^ Vr oscillations. Current Superkamiokande data on single ring events also favors this 
explanation. Data on two ring events can provide a dramatic confirmation of this scenario. 
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FIGURES 

FIG. 1. Plot of Ri vs ^ for S31 = 0.001 eV^ and ^ = (continuous line), = 30° (dotted line) 

and (p = 45° (dashed line). 

FIG. 2. Plot of R2 vs ip for 631 = 0.001 eV^ and ^ = (continuous line), (/) = 30° (dotted line) 
and (f) = 45° (dashed line). 

FIG. 3. Plot of R3 vs tp for ^31 = 0.001 eV^ and ^ = (continuous line), ^ = 30° (dotted line) 
and (f) = 45° (dashed line). 

FIG. 4. Region uKp—ip plane allowed by the Kamiokande constraint on R3. 
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